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The diffuse high-energy gamma-ray emission of the Milky Way arises from interactions of cosmic- 
rays (CRs) with interstellar gas and radiation field in the Galaxy. The neutral hydrogen (H I) 
gas component is by far the most massive and broadly distributed component of the interstellar 
medium. Using the 21-cm emission line from the hyperfine structure transition of atomic hydrogen 
it is possible to determine the column density of H I if the spin temperature (Tg) of the emitting 
gas is known. Studies of diffuse gamma-ray emission have generally relied on the assumption of a 
fixed, constant spin temperature for all H I in the Milky Way. Unfortunately, observations of H I 
in absorption against bright background sources has shown it to vary greatly with location in the 
Milky Way. We will discuss methods for better handling of spin temperatures for Galactic diffuse 
emission modeling using the Fermi-LAT data and direct observation of the spin temperature using 
H I absorption. 



INTRODUCTION 



T using 



The diffuse Galactic emission (DGE) arises from in- 
teractions of cosmic-rays (CRs) with interstellar gas 
and radiation field in the Galaxy. Due to the smooth 
nature of the interstellar radiation field and the CR 
flux after propagation, the fine structure of the DGE 
is determined by the structure of the interstellar gas. 
Getting the distribution of the interstellar gas correct 
is therefore crucial when modeling the DGE. 

It is generally assumed that Galactic CRs are ac- 
celerated in interstellar shocks and then propagate 



N HI (v,T s ) = -log 1- 



2007 



throughout the Galaxy (see e.g. Strong et al 
for a recent review.). In this paper, CR propaga- 
tion and corresponding diffuse e mission is calculated 
using the GALPROP code (see |Strong et al]|2004 
and references within.). We use the so-called conven- 



tional GALPROP model (Strong et al. 2004), where 



the CR injection spectra and the diffusion parame- 
ters are chosen such that the CR flux agrees with the 
locally observed one after propagation. The gas dis- 
tribution is given as Galacto-centric annuli and the 
diffuse emission is calculated for those same annuli. 
The distribution of H I is determined from the 21-cm 
LAB line survey ( |Kalberla et al.|2005[ ) while distribu- 
tion of molecular hydro gen, H 2 , is found usi ng the CO 
(J = 1 — > 0) survey of Dame et al. (2001) assuming 
N R2 = X co (R)Wco- 

While converting observations of the 21-cm H I line 
to column density is theoretically possible, it is not 
practically feasible. To correctly account for the op- 
tical depth of the emitting H I gas, one must know 
its spin temperature, T$ (see e.g. |Kulkarni and Heiles 
1988). Under the assumption of a constant Tg along 



the line of sight, the column density of H I can be 
calculated from the observed brightness temperature 



T 

Ts-Ti 



l>9 



TsC, 



(1) 



where Tb g is the background continuum temperature 
and C = 1.83 x 10 18 cm" 2 K (km/s)' 1 . The assump- 
tion of a constant Tg along the line of sight is known 
to be wrong for many directions in the Galaxy (see e.g. 
Dickey et al.|2009| . The Tg values derived in this pa- 
per are therefore only a global average and should not 
be taken at face value. Figure [T] shows how changing 
Tg affects Nhi in a non-linear way, mainly affecting 
areas with T close to T$ in the Galactic plane. This 
figure was created under the assumption of a fixed 
T$ for the whole Galaxy that is known to be wrong 
but has been used for DGE analysis from the days of 
COS-B ( |Strong|p85| . Note that for equation Q to 
be valid the condition Tg > T + Tb g must hold. When 
generating the gas annuli, this condition is forced by 
clipping the value of T. 

While the assumption of a constant spin tempera- 
ture Ts = 125 K for the whole Galaxy may have been 
sufficient for older instrument, it is no longer accept- 
able for a new generation experiment like Fermi-LAT 
( Atwood et al. 2009 ) . T his has been partially explored 
for the outer Galaxy in Abdo et al. (2010a). In this 



paper we will show a better assumption for Tg can 
be easily found and also show that direct observations 
of Tg using absorption measurement of bright radio 
sources are needed for accurate DGE modeling. 



II. METHOD 

We assume the source distribution of CR nuclei and 
electrons are the same. CR propagation is handled 
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FIG. 1: The ratio 7V H /(125K)/N H i(200K) in Galactic co- 
ordinates. The figure clearly shows the non-linearity of the 
correction that can be as high as a factor of 2 in this case. 



by GALPROP and we use the conventional model 
so that after the propagation the CR spectra agree 
with local observations. The GALPROP diffuse emis- 
sion is output in Galacto-centric annuli, split up into 
different components corresponding to different pro- 
cesses (bremsstrahlung, 7r°-decay, and inverse Comp- 
ton scattering). To allow for radial variations in CR 
intensity we perform a full sky maximum likelihood 
fit, preserving the spectral shape of each component. 
We allow for one global normalization factor for the 
electron to proton ratio. Additionally, we also allow 
for radial variation in the Xqo factor. This accounts 
for uncertainties in the CR source distribution and 
Xco factor. 

The maximum likelihood fits were performed on 



the whole sky using the GaRDiAN package (Abdo 
et al.|[2009 ) after preparing the Fermi-LAT data with 
the science tools. We use the same dataset as lAbdol 



et al. (2010b) that has special cuts to reduce CR 



background contamination compared to the standard 

In addition to 



event selection (Atwood et al. 



the DGE model, we also include all sources from the 



1 year Fermi-LAT source list (Fermi LAT collabora- 
tion|2009 ) and an isotropic component to account for 



2009) 



EGB emission and particle contamination. This fit is 
performed for different assumptions of T$ and a like- 
lihood ratio test is used to compare the quality of the 
fits. 



III. RESULTS 



approximation, apart from the fact that observations 
of the interstellar gas have shown it to be wrong, is 
that the maximum observed brightness temperature 
in the LAB survey is ~150 K which is greater than 
our best fit global Ts- This is solved by clipping the 
observations when generating the gas annuli, which is 
not an optimal solution. A different possibility is to 
use the assumption 



T s = max(T s>1 



+ AT S ). 



(2) 



Here, T max is the maximum observed brightness tem- 
perature for each line of sight. This ensures Ts is 
always greater than T. Scanning the values of Ts iin ; n 
and ATs with a step size of 10 K and 5 K, respectively, 
gives us a maximum likelihood for Ts,min = 110 K and 
ATs — 10 K. While this assumption still does not ac- 
count for the complexity of the interstellar medium, 
the log likelihood ratio between the best fit linear re- 
lation model and the best fit constant Ts model is of 
the order of 1000, a significant change. 

The most accurate Ts estimates come from observa- 
tions of H 1 in absorption against bright radio sources. 
We gathered over 500 lines of sight with observed Ts 
from the literature ( Dickey et al.||2009| |Heiles and 
Troland|[20?53{ |Strasser and Taylor||2004| h This covers 
about 0.2% of the pixels in the LAB survey, allowing 
for accurate column density estimates only in those 
pixels. After taking our best fit linear relation model 
and correcting the pixels with known Ts the fit was 
redone for the whole sky. Note that we did not change 
the values of Ts, m in and AT5. The log likelihood ratio 
of -105 tells us that this model is worse than the best 
fit linear relation. This is not unexpected, since the 
gamma rays are generated from CR interactions with 
the gas and if the gas distribution is wrong, we won't 
get the correct CR distribution from the fit. To limit 
the uncertainty involved with the linear relation as- 
sumption, we did another fit, limiting ourselves to the 
region -10° < \b\ < 10°, 15° < I < 165° that covers 
the observations made in the Canadian Galactic plane 
survey (CGPS) where the density of Ts observations 
is the highest and is large enough to get a good fit 
to the LAT data. The fit in this region results in a 
log likelihood ratio of 28 indicating a statistically sig- 
nificant improvement in the fit. This is despite the 
observed Tg lines of sight only covering 25% of the 
fitted region and the values of Ts !m i n and ATs n °t 
being adjusted after correcting for known Ts values. 



The simplest assumption is that of a constant Ts 
for the whole Galaxy and it deserves some attention 
for historical reasons. It will also serve as a baseline 
model for comparison with other assumptions. To get 
an approximation for the best model, we scan Ts from 
110 K to 150 K in 5 K steps. Our results show that 
Ts = 130 K gives the maximum likelihood for this 
setup. 

One of the problems with the constant global Ts 



IV. DISCUSSION 

Our small exercise here has shown that for accu- 
rate DGE modeling we need to know more about the 
distribution of gas in the Galaxy, especially the H I 
distribution. The standard constant Ts assumption 
is not sufficient for current instruments and small ad- 
justments cause large differences in the quality of the 
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resulting model. We also show that direct observa- 
tions of Ts help in creating a better model of the 
DGE. Unfortunately, direct observations of Tg are dif- 
ficult since they require high resolution telescopes and 
bright radio continuum sources. Some assumptions 
will therefore have to be made for the regions in be- 
tween bright radio sources. 

It must be stated here that all of the above results 
are model dependent. The Fermi- LAT data can only 



provide us with the intensity of gamma-rays from a 
particular direction of the sky. Uncertainties in our 
modeling of contribution other than those directly re- 
lated to the H I distribution will affect the value ob- 
tained for Ts- We are currently studying the system- 
atic effects this will have on our results. We also note 
that even for the best fit models, the residuals show 
signs of structure, strongly indicating our models are 
less than perfect. 
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